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5h . Two-gluon radiative corrections to the h — > cii>i decay width have been computed analytically 

■<^ ' as an expansion in terms of j^^ <C 1 in the kinematical limit of zero lepton invariant mass. The ob- 

^"^ ' tained results match smoothly with a previously known expansion around ( 1 — ?7^ ) ^ 1 • Together 
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they describe the process h — > clv^ for all mass ratios ^^. 



I. INTRODUCTION 

Inclusive semileptonic decays of the 6-quark can be described with a high degree of theo- 
retical precision because nonperturbative effects are suppressed |l| and the 6-quark mass is 
sufficiently large for the perturbative QCD expansion. However, the perturbative calcula- 
tions are challenging because of the presence of massive propagators of b- and c-quarks. In 
particular, the 0{al) corrections to the rate of the decay b — > ciui have not yet been calcu- 
lated, although a rather reliable estimate was obtained in ,2]. The need for such corrections 
has been pointed out in connection with a precise determination of the CKM parameter Vcb 
(see, for example, [3|). Such a calculation requires a generalization of the four-loop studies 
of /i — > eUe^ii [4] and b —>■ uiu£ [5] to the case of a massive charged particle in the final state 
and is thus very difficult. One possible approach is to treat that final-state quark mass as 
a perturbation. In this paper we examine this approach, not with the full problem of the 
total decay rate, but in one special kinematical configuration: the decay b — > ciue with the 
lepton and neutrino escaping with a zero invariant mass, q^ = 0. We treat the mass ratio 
^^ as a small parameter and find several terms of expansion of dr{b -^ ciui)/dq'^ at g^ = 0. 
This work builds on the recent study [£] of the decay t —>■ Wb, where an expansion in the 
W mass was constructed. 

In the limit g^ = 0, massless leptons are produced at rest with respect to each other and 
the corresponding decay rate may be represented as a product of b decay into a fictitious 
real massless VT-boson and a consequent W decay. The decay amplitude is 
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where JT''^ is the quark current. The decay rate oi b -^ Xiue may be represented as follows: 
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Thus dT{b -^ dMi)/dq^ at g^ _, q differs from dV{b -^ cW*) with m{W*) = by only a 
constant factor. In this paper we focus on the decay width r(6 ^ cW*) and integrate over 
final states containing up to two gluons. 



To treat virtual and real gluon emission consistently, we take advantage of the optical 
theorem which connects the 6-quark decay rate to the imaginary part of the 6-quark self- 
energy operator: T{b -^ X) = — Im V{b -^ X -^ b). In this approach the calculation of 
(9(ag) corrections to r(6 -^ cW*) requires the evaluation of 3-loop diagrams instead of the 
4-loop diagrams of the corresponding b —>■ ciug decay. 

Several years ago |2], the b — >• cW* decay rate was found from an expansion around the 



so-called zero recoil limit where ( 1 — ^^^ 

nit 



<^ 1. Our present expansion around ^^^ -C 1 is 



complementary and we will find that, together with [2], we now know the differential rate 
at all values of the c-quark mass. 

This paper is organized as follows. In the next section we introduce the notations, discuss 
the gauge-invariant contributions to the decay rate, and present some technical details of 
the calculation. In Section UTTl we present the results and combine them with those of Ref. 
to cover the full range of final state quark mass. 



II. ©(a^) CORRECTIONS TO 6 ^ cW* DECAY WIDTH. 

The rate oi b —>■ cW* decay may be written as an expansion in the strong coupling 
constant ag whose coefficients are functions of p = ^^: 
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Here Fq = ^ ^ "^*' is the result corresponding to rric = 0. The second-order correction X2 
may be written as a sum of finite, gauge-invariant combinations: 



X2 = TrNlXl + TrNhXh + TrNcXc + CfXa + CaX_ 
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Nl represents the number of massless quarks (3 in this context) while Nh and Nc label the 
contributions of b- and c-quarks, respectively. The top quark contribution is suppressed by 
( ^ j and we neglect it. In S't/(3), the color factors are Tr = |, Cp = |, and Ca = 3. 
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FIG. 1: Examples of diagrams contributing to 0{a1) corrections to 5-decay. 
A. Calculation 

To find Xl, Xh, Xc, Xa, and X^a of Eq. (jH)), we need to consider 39 three-loop diagrams 
such as those in Fig. ^ along with 19 one- and two-loop renormalization contributions. The 
contributing diagrams depend in general on two scales: nib and rric- To account for phe- 
nomena at different scales properly, we apply the well-known asymptotic methods. Table H] 
presents contributing momentum regions in one example of a three-loop double-scale topol- 
ogy. In each region loop momenta are either "hard" (|A;| ^ rric) or "soft" (|/i;| ~ fric), with 
\p\ = 771}, being a hard momentum, allowing us to Taylor expand the propagators and thus 
reduce the number of scales to one. For example. 

Table |l] illustrates the asymptotic expansion process for one of the integrals encountered 
in this work. In this example. Region 4 involves so-called eikonal integrals, featuring the 
propagator 2p(A;3 — ^2) + ^0. Although seemingly double-scale, such integrals only multi- 
plicatively depend on the external momentum p. Care should be taken with eikonal regions, 
since the integral value may depend on the sign of the contour fixing term iO. 

The large number of resulting integrals can be reduced to a small set of "master integrals" 
(see Ref. [7| for an example of a solution algorithm along with references to earlier work). 



Most of the master integrals used in this paper have been described in 



. An additional 



master integral corresponds to the topology in Fig. |2l for completeness we present here the 



TABLE I: Example of the asymptotic expansion of a double-scale topology. In the figures, dotted, 
thin solid, and thick solid lines represent massless, soft-scale massive, and hard-scale massive 
propagators, respectively. 
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(^ai,a2, as, O4, O5, Og, ay, agj — J [x]''i[2]"2[3]a3[4]a4[5]a5[6]''6[7]''7[8]''8 > 

[1] = (fci - Pf, [2] = kl + ml [3] = q + ml [4] = ^ + ml [5] = {p + h - hf + mg, 
[6] = {p + k3- k2f + ml [7] = {k2 - ki)l [8] = (fca - ^2)^, {p^ = -ml) 

.■' .■' '•. '•. Region 1: |ki|, |k2|, jksl » rric (3- loop single-scale topology) 

'■■■..^-~ [2] - kl [3] -. kl [4] - ki 



Region 2: |ki|, jksl » rric, |k2| ~ rric (vacuum bubble x 2-loop topology) 

[2]^kl[4]^kl[6]^ip + ksr + ml[7] 



kU 



Kg 



Region 3: |ki|, |k2| 3> rric, Iks] ~ rrrc (vacuum bubble x 2-loop topology) 

[2] ^ kl [3] ^ kl [8] ^ kl [5] -^{p- k^f + ml [6] ^ {p - k^f + mf 



Region 4: |ki| ^ rric, |k2|, jksl ~ rrrc (eikonal topology x 1-loop topology) 

[2] ^ kl [5] ^ (p - kif + ml [6] ^ 2p(A:3 - ^2)^ + iO, [7] ^ kf 



FIG. 2: The new master integral of Eq. Q. The solid and dashed lines correspond to massive and 
massless propagators, respectively. 



result using the same notations: 
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where J-' = , .^ya is a common loop factor and D = 4 — 2e is the convention used for 
dimensional regularization . 



III. RESULTS 




FIG. 3: Matching of expansions around p = (thick hne) and p = 1 (thin line). The c-quark 
contribution Xc interpolates between heavy (Xh) and light (Xl) quark contributions. 



Our results for the contributions to Eq. (jH)) are obtained as series in p = ^. We have 
obtained terms up to p^ for the expansions and present here contributions to X2 through 
terms of order p^: 
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We have used the MS definition of a^ normahzed at p = ruf,, and the pole mass mi,. 



These expansions may be directly compared to expansions around p 
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as follows: 
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In |2[ corrections from h- and c-quarks were lumped together in A//. Here we divide them 
up by separating the c contribution in A^ (Eq. (jAljl ). As expected, Xc reaches Xl in the 
limit p ^ and Xh in the limit p -^ 1. 

Plots on Fig. El present X^, Xh, Xc, Xa, and Xma calculated to 0(p^°) around p = 0, 
and expansions of the corresponding functions from ^ (Eq. (Al) of that work), according 
to Eqs. (fT^ - (fTH|) . calculated through (9((1 — p)^^). It is sufficient to account for terms up 
to O(p^) to reach the relative accuracy of 1% at the realistic value of rric/mb ~ 1/3. 

For convenience we provide results for a numerical fit providing accuracy better than 0.01 
for Xl, Xc, Xa, and Xjsia, and better than 10~^ for X^ for < p < 1: 



Xl/{1 - pf ^ 2.872 + 6.849p - 17.00p2 + 22.56p3 - 26.92p^ + 13.16p^ 

Xn/il-pf ~ -0.06361 -0.1902p-0.2378p^-0.1733p^-0.09828p^ 

Xc/il - pf ^ 2.882 - 0.9432p - 14.31p2 + 25.00p3 - 18.49p^ + 5.113p^ 

Xa/(1 - pf ^ 3.531 + 1.305v/p + 0.1496p - 13.76p2 + 49.64p^ - 57.77p^ 

+ 33.69p^ - 7.746p^ 

X^a/(1-p)^ ~ -8.090 -1.696v/p- 12. 77p + 20.35p2 + 8.257p^-43.23p^ 

+ 58.52p^ - 29.09p^ 
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IV. CONCLUSION 

Fig. El illustrates existing expansions |3, M, ISi of semileptonic quark decays in various 
kinematic configurations. Analytic expressions are known along the zero recoil line and 
in all corners of the triangle. With the approach demonstrated in this paper and with 
improving computational resources, it is becoming feasible to calculate the complete 0{a'^) 
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FIG. 4: Kinematic boundaries of semileptonic heavy quark decays. The dotted arrow shows the 
expansion presented in this paper. Previously known expansions are indicated with sohd arrows. 
Integration over the dotted line would correspond to integration over the leptonic phase space in 
the decay b -^ civ a. 

correction to the decay rate Yih — > civ) as an expansion around p = 0. Like in the present 
study, the most challenging hard-scale regions will contribute only to even powers of p, with 
odd powers originating from factorized regions. Assuming this series will converge similarly 
to the expansion presented here, a 5% accuracy of the complete correction at realistic value 
P = \ will require calculating terms through p^ of the most difficult diagrams. An extension 
to (^ will require evaluating only factorized diagrams and will likely to improve the accuracy 
to the level of 3%. This is a challenging but feasible task. 

In the future, if the need arises, these expansion techniques could be applied for computing 
precision rates of other heavy colored particle decays, e.g. squarks, accounting for mass- 
dependent effects. 
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APPENDIX A: SOFT QUARK CONTRIBUTIONS 



The contribution of c-quark loops expanded in 5 = 1 — ^^ was calculated as part of the 
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study in Ref. |2] but not explicitly shown. Here, for completeness, we present that result 
through terms 0{S^): 
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